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; ABSTRACT 

The evolution of helium stars with masses of 1.5 - 6.7 Mq in binary systems with 
lO ' a 1.4 Mq neutron-star companion is presented. Such systems are assumed to be the 

remnants of Be/X-ray binaries with B-star masses in the range of 8 - 20 Mq which 
underwent a case B or case C mass transfer and survived the common-envelope and 
spiral-in process. The orbital period is chosen such that the helium star fills its Roche 
lobe before the ignition of carbon in the centre. We distinguish case BA (in which 
mass transfer is initiated during helium core burning) from case BB (onset of Roche- 
^s^' I lobe overflow occurs after helium core burning is terminated, but before the ignition 

. of carbon). We found that the remnants of case BA mass transfer from 1.5 - 2.9 Mq 

' helium stars are heavy CO white dwarfs. This implies that a star initially as massive as 

fSJ I 12 Mq is able to become a white dwarf. CO white dwarfs are also produced from case 

■ BB mass transfer from 1.5 - 1.8 Mq helium stars, while ONe white dwarfs are formed 



> 



^ from 2.1 - 2.5 Mghelium stars. Case BB mass transfer from more-massive helium 

Q^, stars with a neutron-star companion will produce a double neutron-star binary. We 

I ' are able to distinguish the progenitors of type lb supernovae (as the high-mass helium 

O I stars or systems in wide orbits) from those of type Ic supernovae (as the lower-mass 

helium stars or systems in close orbits) . Finally, we derive a " zone of avoidance" in the 
^ ' helium star mass vs. initial orbital period diagram for producing neutron stars from 

helium stars. 
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U ■ Ic 



1 INTRODUCTION 

A binary pulsar with a neutron-star or a heavy white-dwarf 
companion has long been considered to originate from a he- 
lium star in a binary system with a neutron-star companion. 
The latter system is a descendant of a high-mass X-ray bi- 
nary (HMXB), in which the companion of the neutron star 
loses its mass through wind mass loss or a mass-transfer 
phase, exposing its helium core (e.g. Bhattacharya & van 
den Heuvel 1991). Although the existence of such a system 
is only found in Cyg X-3 (van Kerkwijk et al. 1992, 1996), a 
detailed study of the evolution of helium stars in binary sys- 
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tems with a compact companion is important as the systems 
form the bridge between the evolution of X-ray binaries and 
the formation of double compact-object binaries. 

The evolution of a helium star in a binary system has 
been studied e.g. by Savonije, de Kool & van den Heuvel 
(1986) who evolved a 0.6 Mq non-degenerate helium star 
with a 1.3 Mq compact companion. With a short orbital 
period of P = 37™, Roche-lobe overflow (RLOF) takes 
place during helium core burning. Ergma & Fedorova (1990) 
evolved helium stars with masses of 0.5, 0.766, and 1 Mq. As 
the companion star, they took white dwarfs with the same 
range of masses in a combination such that the systems have 
a mass ratio of 0.5 < A/hc/Mwd < 2. The periods are also 
so short, of 26™ 2 - 62™ 6, that the helium stars transfer mass 
to the companion while they are still burning helium in the 
centre. A similar work was also carried out by Tutukov & Fe- 
dorova (1990). A study of systems with more massive helium 
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stars in wider orbits has been carried out by e.g. Delgado & 
Thomas (1981) who considered helium stars with masses of 
2, 2.7, 3.3, and 4 Mq with a massive main-sequence com- 
panion. The helium star fills its Roche lobe after helium is 
exhausted in the core. Habets (1986a) evolved a 2.5 Mq he- 
lium star with a 17 Mq main-sequence companion in a wide 
orbit (P = 20^25) such that mass-transfer phase occurs dur- 
ing first carbon burning convective shell. 

In this work we study the evolution of helium stars with 
masses of 1.5 - 6.7 Mq with a 1.4 Mq neutron-star com- 
panion. Such systems are assumed to be the remnants of 
Be/X-ray binaries with masses in the range of 8 - 20 Mq 
which underwent mass transfer as case B (RLOF is initiated 
during hydrogen shell burning) or case C (during helium 
shell burning). As the result of the large mass ratio, mass 
transfer from the Be star to the neutron star is dynamically 
unstable and the two components are embedded in a com- 
mon envelope (CE), leading to the spiraling-in of the neu- 
tron star in the envelope of the Be star. We assume that the 
system survived the common envelope and spiral-in process. 
We are interested to investigate the final fate of the systems 
(whether they will become white-dwarf/neutron-star or dou- 
ble neutron-star binaries) , the type of supernovae (SN) they 
might produce, as well as to study which systems are stable 
to RLOF. A study on the similar range of mass (2-6 Mq 
helium stars) has been carried out by Avila Reese (1993). 
With a Roche radius of 0.6 and 0.7 Rq, in that study RLOF 
takes place after helium core burning is terminated and the 
calculations were done up to the ignition of carbon in the 
centre. Apart from the slightly wider range of mass in our 
study compared to the latter work and a larger range of or- 
bital periods used in our work, we also are able to follow the 
evolution of the helium stars to more advanced evolutionary 
stage, i.e. beyond the carbon ignition, which enables more 
detailed conclusions on the fate of the systems. We take 
6.7 Mq as the upper limit of our calculation because more 
massive helium stars undergo a dynamically unstable mass 
transfer. 

In Sect. ^ we describe the possibility for the formation 
channel of a helium star in a binary system with a neutron- 
star companion, and some basic assumptions used in the 
calculation of the orbital evolution. The computational code, 
the input parameters used in the code, and the calculation 
that constrains the initial mass of the helium star and the 
initial period are described in Sect. ^ We discuss the results 
in two sections; Sect. ^ for the cases in which RLOF takes 
place during helium core burning and Sect, ^for that during 
helium shell burning. Our conclusions are given in Sect. ^ 



2 THE EVOLUTION OF HELIUM STAR IN 
BINARY SYSTEM 

2.1 Formation of a helium star with a 

neutron-star companion in a binary system 

An HMXB consists of a massive OB-type star and a compact 
object as the X-ray source. There are two classes of HMXB, 
i.e.: (i) OB-supergiant (standard massive) X-ray binaries, in 
which the X-ray source is powered by wind accretion from 
its almost Roche-lobe filling OB-supergiant companion, and 
(ii) Be/X-ray binaries, in which the X-ray source is powered 



by accretion of material when the compact object moves 
through the dense disc of its unevolved Be-star compan- 
ion at periastron passage. The first class is found in close, 
almost circular orbits with companion stars more massive 
than 20 Mq; while the latter in wide, eccentric orbits, with 
8-20 Mq companion stars. 

As the companion starts to fill its Roche lobe, a runaway 
mass transfer to the compact object will follow due to the 
high mass ratio of the system, leading to a CE and spiral-in 
phase. The outcome of this phase is a system consisting of 
the helium core of the companion and a neutron star. The 
helium star evolves and, in turn, might fill its Roche lobe. If 
the helium star is more massive than 2.2 Mq (Habets 1986a) 
and it does not further lose mass, it will collapse and produce 
a neutron star; and the system becomes a double neutron- 
star binary. We refer to Bhattacharya & van den Heuvel 
(1991) for review on the so-called "standard scenario"^ 

An alternative to the standard model was proposed by 
Brown (1995). This new evolutionary scenario suggests that 
a double neutron-star binary is a descendant of a double 
helium-star binary, which is produced from a binary sys- 
tem with components of very similar masses. This twin- 
mass model was proposed because of the argument that a 
neutron star would be able to accrete matter at high rate 
during the CE and spiral-in process and will then become 
a black hole (e.g. Chevalier 1993), which would exclude the 
possibility that double neutron stars would result from the 
CE evolution of HMXB (the standard scenario). Neverthe- 
less, whether or not a neutron star in a CE will indeed go 
to highly super-Eddington accretion is still very uncertain 
(see e.g. Chevalier 1996). Therefore, there is still room for 
the formation of systems consisting of a helium star with a 
neutron-star companion like in the standard scenario. Strong 
accretion onto the neutron star can be avoided, for exam- 
ple, if the donor star is very extended (Chevalier 1993) and 
if rotational effects are taken into account (Chevalier 1996). 



2.2 Orbital evolution 

The orbital angular momentum of a binary system which 
contains a helium star of mass and a neutron star of 

mass Mns is given by 



r 



M 



aG 



(1) 



where M = Mhc + Mns is the total mass of the system, a 
is the orbital separation, and G is the constant of gravity. 
The rate of change in orbital separation of the system is 
expressed by 



Jorb VA^Ho A//nS / M 



(2) 



The total change in orbital angular momentum is assumed 
to be affected by gravitational-wave radiation and by the 



t An evolutionary channel for the formation of double neutron- 
star binary similar to the standard scenario was also proposed 
by Tutukov & Yungelson (1973). They predicted a helium-star 
mass of 2 Mq as the limiting mass for white-dwarf/neutron-star 
formation. 
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loss of mass with angular momentum from the system, i.e. 



J orb 
J orb 



(3) 



where the rate of change of orbital angular momentum duo 
to gravitational-wave radiation is given by (Landau & Lif- 
shitz 1958) 



Jgwr _ 32 MhcMnsM 

Jorb 5 c5 a" 



(4) 



with c the speed of light in vacuum. 

The system is assumed to evolve non-conservatively. 
Mass loss from the system causes loss of orbital angular 
momentum which is taken into account as (van den Heuvel 
1994; Soberman, Phinney & van den Heuvel 1997) 



Jml 

J orb 



1 + q Mhc 



(5) 



where q = A/hc/Mns is the mass ratio; a is the fraction of 
mass lost from the helium star in the form of fast isotropic 
wind, and {3 is the fraction of mass ejected isotropically from 
the vicinity of the neutron star. During the detached phase, 
we took a = 1 and it is assumed that the neutron star does 
not accrete matter from the stellar wind [13 = Q)\ all matter 
released from the system in this case carries the specific or- 
bital angular momentum of the helium star. The stellar wind 
mass loss is approximated by eq. (2) in Wellstein & Langer 
(1999, cf. Hamann, Schonberner & Heber 1982; Hamann, 
Koesterke & Wessolowski 1995) multiplied by a factor of 
0.5 to provide a better fit with the most recent observed 
mass-loss rate of Wolf-Rayet stars (Nugis & Lamers 2000; 
Nelemans & van den Heuvel 2001), i.e. 



Ml 



He ,wind 



\1.5 



2.81O-"(L/L0)^ 
4.O1O-3^(L/L0)«'«, 



logL/L0 > 4.5 
log L/Lq < 4.5 



(6) 



This means our adopted mass loss rate in the upper lumi- 
nosity range is one quarter of the Hamann et al. (1995) rate. 
During RLOF, the helium star transfers mass with a mass- 
transfer rate which is defined as 



Ml 



He, RLOF ~ 



-10 X max 



0, In 



Rl 



(7) 



where _Rl is the helium star's Roche radius given by Eggle- 
ton (1983) as 



0.49 a 



0.6 + q-2/3 ln(l + gi/3) 



(8) 



During the mass-transfer phase (even on a nuclear timescale) 
the wind mass loss contributes at most 5 per cent of the 
mass-loss rate. It is assumed that the neutron star accretes 
matter up to its Eddington limit for helium accretion (3 x 
10~* M0yr~^); the rest of the transferred matter is lost 
from the system with the specific orbital angular momentum 
of the neutron star. 

In the limiting case when the mass-transfer rate owing 
to RLOF is much higher than both the maximum accretion 
rate and the wind mass-loss rate (M > 10~® MQyr"^), we 
have p 1 and Jgwr can be neglected. This is generally the 
case during case BAB and case BB mass transfer discussed 
below. Then eq. (^ leads to: 



P 
P 



3g2 - 2g - 3 Mne 



(9) 



1 + q Mhs 

From this it can be seen that P has a minimum when 
q — 1.39, i.e. Mne ~ 1.95. For q > 1.39, P initially decreases 
as a result of RLOF, and increases again after q < 1.39. This 
behaviour can be seen in Sect. I and H 



3 THE METHOD OF CALCULATION 

3.1 A brief description on the numerical code 

We used the Eggleton stellar-evolution code (Eggleton 1971, 
1972, 1973). This code uses a self-adaptive, non-Lagrangian 
mesh-spacing which is a function of local pressure, tempera- 
ture, Lagrangian mass and radius. It treats both convective 
and semi-convective mixing as a diffusion process and finds a 
simultaneous and implicit solution of both the stellar struc- 
ture equations and the diffusion equations for the chemical 
composition. The most recent update, including the nuclear 
reaction network, is described in Pols et al. (1995). 



We evolved zero-age main-sequence stars (in Sect. 3.2) 
assuming a chemical composition of {X — 0.70, Z = 0.02) 
and using a mixing- length parameter of a = I /Hp = 2.0. 
Convective overshooting is taken into account in the same 
way as in Pols et al. (1998) using an overshooting constant 
5ov = 0.12. 

The helium stars were evolved assuming a chemical 
composition of {Y = 0.98, Z = 0.02). All calculations were 
done without enhanced mixing (STD models in Pols 2002). 
We refer to the latter paper for the evolution of single helium 
stars related to our calculations. 



3.2 Initial conditions 

3.2.1 Progenitor systems of helium star and neutron star 
binaries 

To obtain estimates of the mass of the helium star and the 
initial period, we evolved massive main-sequence stars with 
masses of 8 - 20 Mq. If the radius of the star is assumed 
to be equivalent to its Roche radius when RLOF is initiated 
{R ~ Pl), then the separation of the binary system at the 
onset of the mass-transfer phase, ai, for a given radius R can 
be calculated with eq. (^). With the energy equation for CE 
evolution (Webbink 1984; de Kool 1990) 

Of McoroM2 1 

fli 



A/do 



Ma + 2Menv/(?7CEArL) 



(10) 



we calculated the final separation, Of. Here Mdonor ~ Afcore + 
Monv, Mcore, and AJenv are the masses of the donor star, its 
helium core and envelope, respectively. Due to the stellar 
wind mass loss (de Jager, Nieuwenhuijzen & van der Hucht 
1988; Nieuwenhuijzen & de Jager 1990), .Afdonor might be 
smaller than Mzams, the mass of the star on the zero-age 
main sequence (ZAMS). During the red supergiant phase, 
the stellar mass is reduced by 3 - 19 per cent (in a 8 Mq 
and a 20 Mq stars, respectively). A/2 is the mass of the 
companion, which we assumed to be a 1.4 Mq neutron star. 
rL = Rh/ai is the dimensionless Roche radius of the donor 
star. We took ??ce, the so-called CE efflciency parameter. 
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Figure 1. Mass of the helium core and orbital period following 
a CE evolution resulting from case B and C mass transfers from 
8-20 M0 main-sequence stars. The numbers below the figure 
indicate the corresponding initial ZAMS mass. The full lines sep- 
arate the regions of case BA, BB, and BC mass transfers, taken 
from Pols' (2002) single helium-star calculations. The dotted lines 
mark the maximum period allowed resulting from case B (lower- 
and upper-line, derived from the two extremes of A- value); and 
the dashed line resulting from case C mass transfer. See text for 
an explanation. 



to be one. A, the parameter of the binding energy of the 
envelope to the core, is approximated as in Dewi & Tauris 
(2000). We assumed that the binary will survive the CE 
and spiral-in phase if the final separation is larger than the 
Roche radius of the helium core, i.e. the helium star does 
not immediately fill its Roche lobe. Following Taam, King 
& Ritter's (2000) argument, we assume that the neutron 
star does not accrete a significant amount of matter while it 
is embedded in the donor's envelope. 

The results of the calculations are summarized in Fig. ^ 
where we have plotted the maximum periods (which are 
related to the maximum radius) after case B and case C 
spiral-in. These values are very uncertain, especially for case 
B, because the parameter A used for calculating the post- 
CE separation depends quite sensitively on the choice of the 
mass of the helium core. As argued by Tauris & Dewi (2001) 
approximate lower limits to A and to the core mass are ob- 
tained by assuming that the core mass is the mass where the 
abundance of hydrogen is below 0.1. Considering the bottom 
of convective envelope as the core, results in approximate 
upper limits of A and the core mass. The latter criterion 
corresponds roughly to the mass where Xh = 0.3, for the 
maximum radius during case B from AfzAMS < 13 Mq. The 
dotted lines in Fig. |l| mark the maximum periods for case 
B evolution, calculated with the two extreme A-values. For 
case C evolution, only the maximum period (dashed-line) 
corresponding to the lower limit of A is plotted. In stars with 
AfzAMS ^ 14 Mq, helium is ignited before a deep convec- 
tive envelope has developed, and the radius keeps expanding 
during core He burning. At the maximum radius for case B 
(which in this case is not well-defined), the envelope is more 
strongly bound and the maximum period after case B is cor- 
respondingly smaller (see Fig. ^ . In this case, the bottom of 



Table 1. Terminology of the mass-transfer phases 





Case B 


Case C 


Occurence 


H-shell burning 


He-shell burning 


Outcome 


He star 


He star + CO-core 


Helium star evolution 






RLOF during 






• He-core burning 


Case BA 




successive He-shell burning 


Case BAB 




• He-shell burning 


Case BB 


Case CB 


• C-burning 


Case BC 


Case CC 



the convective envelope can no longer be realistically taken 
as the core mass, so the upper dotted line in Fig. hi stops 
at 13 Mq. Because of these uncertainties, to relax the as- 
sumption on the choice of core mass as well as on 77CE (i.e. 
the possibility of applying rycE other than unity), we did the 
calculations for a wide range of periods. 



3.2.2 Studied types of helium star binary evolution 

As Fig. ^ shows, the progenitors of helium star and neu- 
tron star binaries were main-sequence star and neutron star 
systems that evolved according to case B or case C mass 
transfer. The outcome of CE evolution in case B is an un- 
evolved helium star and neutron star, while that in case C is 
an evolved helium star (which has developed a carbon core) 
and neutron star. The systems resulting from case B can be 
so close that they start RLOF during helium core burning. 
This is so-called case BA evolution. Wider systems will start 
RLOF only during helium shell burning around the carbon 
core. This is so-called case BB evolution. The terminology 
of the types of evolution used in this paper is summarized 
in Table |l[ 

Case C mass transfer will result in evolved helium stars 
which, after CE evolution, can only start transferring mass 
after helium core burning is terminated. We are interested in 
RLOF during helium shell burning, which can be called case 
CB evolution. Basically, case BB and CB evolutions result 
in the same type of RLOF from a helium star of about the 
same mass and period, except in A^hc > 4 Mq where the 
wind mass loss plays a major role in reducing the mass and 
changing the structure during helium core burning. We do 
not separately study case CB evolution in this paper, but 
assume that it is similar to case BB evolution. 

A single helium star expands and reaches the maxi- 
mum radius during convective carbon shell burning (Habets 
1986a; Pols 2002). Although they are not less important, we 
did not attempt to follow evolutions in which RLOF is initi- 
ated during carbon core burning or beyond (case BC or CC 
evolution). The overall results for case BA mass transfer are 
shown in Table gI and those for case BB in Tables p and n. 
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Table 2. The binary parameters in case BA mass transfer from 1.5 - 2.9 Mq helium stars: the initial mass and period; the duration, 
amount of mass transferred from the helium star, and the maximum mass-loss rate during the case BA mass transfer; the duration of 
the detached phase; the duration, amount of mass transferred from the helium star, and the maximum mass-loss rate during the case 
BAB mass transfer; the final mass and period, the final mass of helium in the envelope and the final mass of the neutron star. Masses 
are in solar unit, mass-loss rates in Mq yr^^, times in yr and periods in days. * denotes the slow phase of mass transfer. 
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Pc 




^NS.f 






^*RLOF 


AM 










AM 


-^^max 












1.5 


0.050 


4.51 X 10*^ 


0.290 


1.4 X 10- 


. 7 


1.17 X 10° 


8.27 X 10^ 


0.425 


1.4 X 10- 


-6 


0.785 


0.11 


0.0350 


1.544 




0.060 


9.69 X 10^ 


0.042 


8.3 X 10- 


-8 


9.88 X 10^ 


4.96 X 10^ 


0.580 


3.1 X 10- 


-6 


0.864 


0.12 


0.0262 


1.439 


1.8 


0.050 


3.84 X 10'^ 


0.440 


3.3 X 10" 


-7 


9.94 X 10^ 


6.00 X 10^ 


0.513 


2.2 X 10- 


-6 


0.836 


0.12 


0.0276 


1.523 




0.060 


1.08 X 10*^ 


0.138 


2.7 X 10- 


-7 


8.76 X 10^ 


3.48 X 10^ 


0.715 


5.2 X 10- 


-6 


0.933 


0.13 


0.0168 


1.441 


2.1 


0.050 


3.74 X 10'^ 


0.683 


1.2 X 10- 


-6 


9.81 X 10^ 


5.41 X 10^ 


0.556 


2.7 X 10- 


-6 


0.857 


0.12 


0.0277 


1.521 




0.060 


1.16 X 10^ 


0.360 


9.7 X 10- 


-7 


8.07 X 10^ 


3.03 X 10^ 


0.768 


6.1 X 10- 


-6 


0.970 


0.13 


0.0133 


1.443 


2.4 


0.060 


1.37 X 10^ 


0.681 


5.0 X 10- 


-6 


8.13 X 10^ 


3.15 X 10^ 


0.731 


5.4 X 10- 


-6 


0.981 


0.12 


0.0119 


1.449 










1.1 X 10- 
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0.065 


5.16 X 10^ 


0.432 


4.6 X 10- 


-6 


6.25 X 10^ 


1.92 X 10"^ 


0.877 


9.8 X 10- 


-6 


1.089 


0.11 


0.0147 


1.421 










8.6 X 10- 
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2.5 


0.060 


1.48 X 10^* 


0.806 


9.5 X 10- 


-6 


8.29 X 10^ 


3.32 X 10^ 


0.707 


5.1 X 10- 


-6 


0.976 


0.12 


0.0122 


1.453 
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0.065 


6.58 X 10^ 


0.585 


9.8 X 10- 


-6 


6.17 X 10^ 


2.09 X 10^ 


0.821 


8.1 X 10- 


-6 


1.079 


0.11 


0.0109 


1.425 



7.2 X lO-'^ 



2.8 


0.060 


5.50 X 10* 


1.061 


4.9 X 10- 


-5 


1.45 X 10^ 


3.47 X 10^ 


0.713 


5.4 X 10- 


-6 


0.968 


0.12 


0.0091 


1.429 






6.16 X 10^ 


0.049 


1.5 X 10- 


-7 


7.90 X 105 




















0.065 


5.51 X 10-* 


1.032 


5.1 X 10- 


-5 


1.74 X 10*5 


2.81 X 10^ 


0.721 


5.3 X 10- 


-6 


1.041 


0.11 


0.0057 


1.410 


2.9 


0.060 


5.56 X 10* 


1.394 


7.3 X 10- 


-5 


3.92 X 10^ 


4.38 X 10^ 


0.577 


2.4 X 10- 


-6 


0.922 


0.13 


0.0047 


1.415 




0.065 


5.43 X 10* 


1.330 


7.6 X 10- 


-5 


2.29 X 10*5 


3.49 X 10^ 


0.566 


2.0 X 10- 


-6 


0.999 


0.12 


0.0071 


1.412 



4 RESULTS: CASE BA MASS TRANSFER 

Based on the stability of RLOF during helium core burning, 
we devide case BA mass transfer into three groups, i.e. (i) 
1.5 - 2.1 Mq helium stars, in which RLOF is stable, (ii) 
2.4 - 2.5 Mq helium stars, in which RLOF is initiated by a 
rapid, unstable phase and followed by a slow, stable phase of 
mass transfer, and (iii) 2.8 - 2.9 Mq helium stars, in which 
the unstable, thermal-timescale RLOF causes the temporary 
disappearance of the convective core. 

In general, case BA mass transfer from helium stars 
in our calculations resembles case A evolution from main- 
sequence stars (e.g. Pols 1994; Wellstein, Langer & Braun 
2001). In all helium stars of 1.5 - 2.9 Mq which undergo 
case BA mass transfer, the whole envelope is removed from 
the star, and the remnant, a degenerate CO white dwarf, is 
never massive enough to ignite carbon. 

4.1 Roche-lobe overflow from 1.5 — 2.1 Mq 
helium stars 

During the detached phase the orbit slowly shrinks due to 
gravitational-wave radiation and the stellar radius slowly in- 
creases, until (stable) RLOF starts (point A in Fig. Dur- 
ing the first phase of stable, nuclear-timescale RLOF, the 
star contracts and the energy-production rate in the centre 
decreases as a result of the loss of mass, which causes the 
luminosity and effective temperature to decrease. The star 
remains in thermal equilibrium. The mass of the convective 



core does not grow as in the case of a single star, but de- 
creases slowly until the core disappears when the helium is 
exhausted (see Fig. Most of the transferred matter is lost 
from the system at low specific orbital angular momentum 
which, although counteracted by gravitational-wave radia- 
tion, causes the orbit to slightly expand. Roughly at the 
same time when the abundance of carbon in the centre is 
maximum, the period reaches a maximum value, while L 
and TefF reach their minimum values (the reddest point in 
the HR diagram, point B). After this point, the orbit gradu- 
ally shrinks. We find that the earlier RLOF is initiated (the 
smaller the initial period), the more mass is removed from 
the helium star, the more significant the change in period, 
and the longer the duration of RLOF. 

After the turning point, the star moves to the left in 
the HR diagram, only an insignificant amount of mass is 
transferred from the helium star - this is also the case if the 
orbital period is such that mass transfer starts after the car- 
bon abundance in the centre reaches its maximum. In the 
latter case, there is no maximum period; the period contin- 
uously decreases (while the star shrinks) until the system 
is detached. The maximum mass-loss rate, the duration of 
RLOF, and the amount of mass transferred to the neutron 
star all decrease with decreasing initial mass ratio and with 
increasing initial period. 

As the helium abundance in the core drops below about 
0.1, the star contracts inside its Roche lobe, and the system 
becomes detached (point C). During this detached phase, 
the orbit is again tightened due to gravitational-wave radi- 
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8.6 



8.4 




log Pc (g/om ) 



4.5 




Figure 2. The evolutionary tracks in the HR, diagram (bottom) 
and central density-central temperature plane (top) for case BA 
mass transfer from a 1.5 Mq (dotted-) and a 2.8 Mq helium stars 
(solid-line). See text for an explanation of the labels. 



ation, but the net change in orbit is insignificant compared 
to the initial period. The duration of helium core burning 
is prolonged due to the mass transfer. However, we found 
that if the initiation of RLOF takes place close to the end 
of the helium main sequence, helium core burning proceeds 
marginally faster (by less than 1 per cent) than in a single 
helium star of the same mass. 

After helium is exhausted in the centre and the con- 
vective core disappears, helium burning moves to a radia- 
tive shell around the CO core. The CO core grows in mass 
but decreases in radius, which causes the star to expand 
and move up in the HR diagram. A second phase of sta- 
ble RLOF, which we call case BAB, then takes place (point 
D). The central temperature increases. The mass-loss rate, 
which is higher than that during the helium core burning, 
increases both with initial mass ratio and with initial period. 
The orbit expands in all cases due to the effects described in 
Sect. 



2.2 



At the end of this mass-transfer phase, when al- 
most the entire envelope is removed from the star, the degen- 
eracy border is crossed. The central temperature decreases 
as well as the energy-production rate. Although the core is 
constant in mass, the radius decreases and the star contracts 
inside its Roche lobe. The system is again detached (point 
E). The duration of the second RLOF phase decreases with 
increasing initial period and initial mass ratio. 



0.1 

^ 0.08 
a. 

0.06 



1.5 Mq He + 1.4 Mq NS, P^^^ = 0.05 d 




4x10° 6x10° 8> 

time (yr) 

Figure 3. Case BA. The evolution of a 1.5 Mq helium star 
with a 1.4 M0 neutron-star companion in a very close orbit. The 
upper panel shows the orbital evolution. The second panel shows 
the central abundances. Solid-, dashed-, dash-dotted-, and dotted- 
lines are the helium, carbon, oxygen, and neon abundances, re- 
spectively. The third panel gives the mass-loss rate in MQyr"^; 
wind mass-loss rate during the detached phase, or mass-transfer 
rate during the Roche-lobe overflow phase. The fourth panel 
presents the stellar (solid-) and Roche radii (dotted- line), in solar 
radius. The fifth panel shows the evolution of the stellar interior. 
The solid line represents the total mass of the donor star; the dot- 
ted line the CO core mass (defined as the mass where the helium 
abundance is below 0.1). The dark- and light- shaded areas mark 
the helium-burning convective and semiconvective regions. 



4.2 Roche-lobe overflow from 2.4 - 2.5 Mq 
helium stars 

RLOF is initially unstable and results in a rapid, thermal- 
timescale phase of mass transfer. The rapid loss of mass 
causes a rapid contraction of the star. The mass of the con- 
vective core decreases rapidly when the star adjusts to the 
new total mass. The luminosity and effective temperature 
decrease, as well as the central temperature and the orbital 
period, until the mass-loss rate reaches a maximum value. 
The mass of the convective core then decreases slowly and 
the star contracts slowly. Mass transfer continues stably on 
a nuclear timescale. In a 2.5 Mq helium star, the transi- 
tion between the rapid, thermal-timescale phase and the 
slow, stable mass-transfer phase is marked by a minimum 
in mass-loss rate (see Fig. ^). The mass of the convective 
core reaches a minimum value before it increases and then 
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2.5 Mq He + 1.4 Mg NS, P^^^ = 0.06 d 
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Figure 4. Case BA. Same as Fig. g for a 2.5 Mq tielium star. 
Ttie peak and a dip in the mass-loss rate mark the rapid, unstable 
mass-transfer phase, and the transition to a slow, stable phase, 
respectively. 



0.12 
„ 0.1 
I 0.08 

Q. 

0.06 
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time (yr) 



Figure 5. Case BA. Same as Fig. ^ for a 2.8 Mq helium star. 
Note the constancy in central abundance, a dip in the period, and 
the disappearance of the convective core during the first, unstable 
mass-transfer phase. 



decreases slowly. The further evolution is v ery similar to that 
of the lower-mass stars described in Sect. |4.l| . 

As in the case of 1.5 - 2.1 Mq helium stars, the star 
shrinks inside its Roche lobe before helium is depleted in the 
centre. The star again fills its Roche lobe when the radius 
increases after helium is exhausted in the centre (case BAB). 

4.3 Roche-lobe overflow from 2.8 - 2.9 Mq 
helium stars 

The mass transfer from 2.8 - 2.9 Mq helium stars is also ini- 
tiated (point 1 in Fig. |2|) by a thermal-timescale phase. This 
time, the rapid mass transfer with q ^ 2 causes the star 
to lose about 40 per cent of its mass rapidly. The central 
temperature decreases below the required temperature for 
helium burning so that the energy production in the centre 
ceases (point 2, see also Fig. This, in turn, causes rapid 
contraction of the star and disappearance of the convective 
helium core. The period decreases rapidly. After the maxi- 
mum mass-loss rate is reached, the star contracts inside its 
Roche lobe, the orbit widens due to the loss of mass from 
the system, and the system is detached (point 3). 

After thermal stability is restored (point 4), the convec- 
tive core reappears, helium is reignited, and the star slowly 
expands. Wind mass loss is now that of a ~ 1.7 Mq helium 



star. The orbital period again decreases due to gravitational- 
wave radiation. In a 2.8 Mq helium star in a very close orbit, 
at the moment when the central carbon abundance reaches 
a maximum, the stellar radius reaches its Roche radius, and 
a brief phase of stable RLOF takes place (point 5 to 6). Bi- 
naries of the same mass in a slightly wider orbit (as well 
as those with Mhb = 2.9 Mq) do not undergo this sec- 
ond, slow phase of mass transfer since the maximum radius 
(which relates to the maximum energy-production rate and 
to the maximum central carbon abundance) is still within 
its Roche radius. However, the net effect is the same since 
during the slow phase only an insignificant amount (2 per 
cent) of mass is removed from the helium star. As in the 
case of lower-mass stars, the 2.8 - 2.9 Mq helium stars also 
undergo case BAB mass transfer (point 7 to 8) with a lower 
mass- loss rate than that during the first, rapid RLOF phase. 

In more massive helium stars (Mhc > 2.9 Mq), case 
BA causes the loss of more than half of the stellar mass. 
The rapid mass loss causes the disappearance of the convec- 
tive helium core. However, the star is unable to recover its 
thermal equilibrium. RLOF becomes dynamically unstable 
and the binary will probably coalesce. 
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4.4 The role of gravitational-wave radiation 

Before the star fills its Roche lobe, the orbit is tightened due 
to gravitational-wave radiation. In order to see the effect 
of gravitational-wave radiation, we also studied the case in 
which gravitational-wave emission is turned off. The exclu- 
sion of loss of angular momentum due to gravitational-wave 
radiation will only delay and shorten the first RLOF phase. 
In this case, during the slow stable RLOF phase, the orbit 
expands due to the loss of mass from the system, and there 
is no maximum period at the moment the central carbon 
abundance reaches its maximum. 

The exclusion of gravitational- wave emission in the case 
of higher-mass helium stars (A/ho > 2.9 Mq) with very close 
orbits allows the wind mass loss to govern the loss of angular 
momentum. The orbit is widened such that mass transfer 
will only take place after the exhaustion of helium in the 
centre (case BB). 



5 RESULTS: CASE BB MASS TRANSFER 

We devide case BB mass transfer into two mass groups, i.e. 
(i) 1.5 - 2.1 Mq helium stars, in which RLOF is terminated 
when the degeneracy border is crossed and a heavy, degen- 
erate CO or ONe white dwarf is produced, and (ii) 2.4 - 
6.6 Mq helium stars, in which an ONe white dwarf or a 
neutron star is the outcome. 



5.1 Roche-lobe overflow from 1.5 
helium stars 



2.1 M 







Case BB mass transfer from 1.5 - 2.1 Mq helium stars re- 
sembles case BAB in the same range of mass. The exhaus- 
tion of helium in the core causes the star to expand and fill 
its Roche lobe. After the degeneracy boundary is crossed, 
it contracts rapidly inside its Roche lobe, and the system 
is detached. The mass-loss rate increases while the duration 
of RLOF decreases with increasing initial period and initial 
mass ratio. 

After the core becomes degenerate but before the sys- 
tem detaches, carbon burning occurs in the radiative shell 
with maximum temperature. In 2.1 Mq helium stars, off- 
centre C-ignition takes place in the shell with maximum 
temperature. At first, carbon burns in a broad radiative shell 
(and marginally in the centre). When the off-centre convec- 
tive shell appears, the whole core expands and the region 
interior to the convective shell cools. C-burning stops com- 
pletely in the inner core, while convective off-centre burning 
continues. The calculations break down shortly after this 
point. It is expected (see Pols 2002 for a discussion) that 
eventually C-burning reaches the centre and a degenerate 
ONe core is formed. In the 1.5 and 1.8 Mq helium stars, 
the temperature is never high enough to cause significant 
carbon burning. 

Before RLOF occurs, the gravitational-wave radiation 
brings the components closer together. In the systems with 
wider orbits (P > 0^*5), the effect of gravitational-wave ra- 
diation is weak so that the period remains constant. The 
evolution of the period is as expected from Sect. 



2.2 



Table 3. The binary parameters in case BB mass transfer from 
1.5 - 2.1 Mq lielium stars: the initial mass and period; the du- 
ration, amount of mass removed from the helium star, and the 
maximum mass-loss rate during the Roche-lobe overflow; the fi- 
nal mass and period, the final mass of the helium in the envelope 
and the final mass of the neutron star. Masses are in solar mass, 
mass- loss rate in Mq yr"-*^, time in yr, and period in days. 



Mi Pi AfRLOF AM 



Mn: 



Mf Pf Mhc,c A^i 



NS.f 



1.5 0.08 3.81 X 10^ 
0.50 1.16 X 10^ 
1.00 8.77 X 10* 
1.8 0.08 2.41 X 10^ 
0.50 6.53 X 10* 
1.00 4.94 X 10* 
0.08 1.51 X 10^ 
0.50 3.22 X 10* 
1.00 2.03 X 10* 



2.1 



0.581 2.2 
0.559 9.8 
0.527 1.5 
0.793 5.6 
0.749 2.4 
0.748 3.8 
0.966 1.2 
0.874 4.8 
0.852 7.6 



X 10" 
X 10" 
X 10" 
X 10" 
X 10" 
X 10" 
X 10" 
X 10" 
X 10" 



" 0.897 0.16 0.0154 1.411 

•5 0.940 0.97 0.0168 1.404 

5 0.946 1.90 0.0164 1.403 

« 0.996 0.16 0.0098 1.408 

5 1.043 0.91 0.0070 1.402 

5 1.046 1.79 0.0067 1.402 

5 1.126 0.14 0.0220 1.404 

^ 1.224 0.73 0.0567 1.401 

5 1.233 1.43 0.0632 1.401 



that RLOF is about to end. Almost the entire envelope is 
removed from the star during RLOF, supporting Delgado 
& Thomas' (1981) resuh. The remnants of 1.5 and 1.8 Mq 
helium stars are degenerate CO white dwarfs with masses of 
0.9 - 1.05 Mq, which are more massive than the remnants 
of case BAB mass transfer from helium stars in the same 
range of mass. This trend resembles Pols' (1994) result for 
case AB and B mass transfers. Helium stars with initial mass 
2.1 Mq produce ONe white dwarfs. 



5.2 Roche-lobe overflow from 2.4 - 6.6 Mq 
helium stars 

The following classification is based on the time of occur- 
rence of the maximum mass-transfer rate; whether it occurs 
far before, just before, or during carbon core burning. This 
can be easily recognized from their tracks in the HR dia- 
gram, as are shown in Fig. ^ The close orbit systems are 
characterized by a small hump (indicated by a star) which 
marks the transition from the rapid to the slow mass-transfer 
phase. The wide orbit systems are characterized by the u 
and down movement (see the right-hand side inset in Fig, 
which is caused by the decreasing and increasing of the mass- 
loss rate due to the appearance and disappearance of the C- 
burning convective shells (except in a 2.4 Mq helium star) . 
The intermediate orbit systems show none of the two char- 
acteristics. 

Before RLOF is initiated, gravitational-wave radiation 
only plays a role in shrinking the orbit in helium stars with 
masses of 2.4 - 3.4 Mq in close orbits during the early phase 
of their evolution. In stars of the same masses in wider or- 
bits, gravitational-wave emission is negligible so that the 
orbit remains constant. Wind mass loss dominates the loss 
of angular momentum and widens the orbit during the gi- 
ant phase of the above range of mass, and during the whole 
evolution of Mhg > 3.7 Mq stars. 



We are not always able to follow the evolution until the 
detachment. However, the drop in mass-loss rate suggests 



© RAS, MNRAS 000, §-0 



The evolution of helium stars in close binary systems 9 



Table 4. Same as Table ^for 2.4 - 6.6 Mq helium stars.* and ** denote the second and third phases of Roche-lobe overflow, respectively. 
The number in the Note column indicates the carbon shell burning in which the calculation stops, a implies the appearance of the shell, 
b the shell burning, and d the disappearance of the shell; c signifies that the calculation goes up to the penetration of the helium-burning 
convective shell to the surface. For 2.4 and 2.5 Mq stars, cb denotes the C-burning and cd the depletion of carbon in the off-centre shell. 
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2.4 


0.08 
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Figure 6. The HR diagram for a 2.8 Mq helium star with a 1.4 Mq neutron-star companion in a close (Pi = 0'?08), intermediate 
(P; = O^S), and wide orbit (Pj = l'*). The triangles mark the onset of (first) mass transfer phase, and the rectangles the (first) maximum 
mass-loss rate. See text for an explanation about the labels. 



5.2.1 Close-orbit systems 

In the close-orbit systems, mass transfer takes place initially 
on the donor star's thermal timescale. The star contracts 
and the luminosity and efTective temperature - as well as 
the period - decrease, until the mass-loss rate reaches a 
maximum, Mmax,!- Toff then increases until the point which 
marks the transition between the rapid- and slow-phase of 
mass transfer (the hump marked by a star in Fig. ^). L 
reaches another minimum, while P also decreases to a min- 
imum (point A). The star then moves to the upper left in 
the HR diagram. 

The appearance of a C-burning convective core (or off- 
centre shell in A/hc < 2.5 Mq stars) causes an increase in 
mass-loss rate to another maximum, Mmax, 2, which varies 
from 33 (in a 5.4 Mq helium star) to 58 per cent (in a 
2.4 Mq helium star) of Mmax.i. After the convective core 
stops growing in mass in Mue < 4 Mq stars, during carbon 
core burning, the mass-loss rate drops, the star moves to 
the left in the HR diagram and shrinks inside its Roche lobe 
(point B, see the left-hand side inset in Fig. ^). The period 
is relatively constant during this detached phase. 



When the carbon abundance in the convective core 
drops below 0.1, the star expands again, moves to the right 
in HR diagram, and the orbit expands .|^ The second phase of 
RLOF occurs after the C-burning convective core has van- 
ished (point C). The mass- loss rate increases rapidly and 
reaches a maximum, Mmax, 3 (point D), when the first C- 
burning convective shell appears. When this shell disap- 
pears, the mass-loss rate drops again. 

In Mhb < 2.9 Mq stars, the radius decreases and the 
components are again detached (point E) . The star moves to 
the left again at constant period until it reaches a minimum 
radius when the second C-burning convective shell appears. 
The star fills its Roche lobe for the third time (point F) 
after the second convective shell has vanished. M increases 
and reaches a maximum value, Mmax, 4- The second RLOF 
lasts shorter than the first one, and the third one shorter 



3 We could not follow the subsequent evolution in A/jje < 2.5 Mq 
stars. However, we expect that these stars will undergo the same 
phenomenon as is found in the more massive ones. 
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than the second, because the evolution of the interior region 
proceeds faster (Habets 1986a). 

In 3.1 - 4 Mq helium stars, RLOF takes place in two 
episodes. After the first C-burning convective shell disap- 
pears, the mass-loss rate decreases to a minimum, but RLOF 
continues without a detached phase. In 3.6 - 4 Mq helium 
stars, the first shrinkage of the star inside its Roche lobe 
takes place later than in the lower-mass stars, i.e. when car- 
bon is being depleted in the core. The radius reaches its 
minimum value after the core has vanished. With the ap- 
pearance of the first C-burning convective shell, the star 
fills its Roche lobe for the second time. 

In Mhb > 4 Mq stars, M also decreases after reaching 
Mniax,2, until carbou is exhausted in the core. The behaviour 
of M subsequent to carbon core burning is erratic, but it 
never exceeds Mmax,i. In this range of mass, the star never 
shrinks inside its Roche lobe. 

Why detached phases occur in one case and not in 
the other, and why the detachments take place in different 
phases of carbon burning for different cases, can be under- 
stood as follows. The start of each convective carbon burning 
phase (both in the core and in subsequent shells) causes the 
core to expand. As a result, the envelope tends to contract 
and the star tends to shrink inside its Roche lobe. In between 
the convective phases, the core contracts and the envelope 
expands again. However, the reaction of the stellar surface is 
delayed by the thermal diffusion timescale between the core 
and the surface. This timescale is roughly constant for each 
star, and does not vary strongly with mass. On the other 
hand, the timescale for C-burning decreases strongly with 
each successive convective phase, as well as with increasing 
mass, and eventually become much shorter than the thermal 
timescale. As a result, the detachments are delayed or even 
absent for advanced phases and/or large masses. 

The transition from rapid to slow mass transfer is not 
found in Mho > 5.4 Mq stars. Stars in this mass range with 
periods shorter than those in Table ^ go through dynami- 
cally unstable mass transfer. 



5.2.2 Intermediate- orbit systems 

In the intermediate-orbit systems, RLOF occurs on the ther- 
mal timescale, the orbit shrinks, and the maximum mass- 
transfer rate, Mniax.i, is reached not long before the appear- 
ance of the C-burning convective core (or off-centre convec- 
tive shell in 2.4 - 2.5 Mq helium stars). The star moves 
afterwards to the upper left in the HR diagram (see Fig. |^. 
Unlike in the close-orbit systems, the mass-transfer rate does 
not increase when the convective core appears. Along with 
the increase in L and TcS, M decreases. 

In 2.4 - 2.9 Mq helium stars, AI decreases until the 
star shrinks inside its Roche lobe during carbon core (or 
off-centre shell) burning, and the system is detached (point 
a, see the middle inset of Fig. ^ and Fig. The period 
increases shortly after Afniax.i is reached, and is relatively 
constant during the detached phase. When the carbon abun- 
dance in the core (or the off-centre shell) drops below 0.1, 
the star reaches the bluest point (In 2.4 - 2.5 Mq helium 
stars, at this point the off-centre shell penetrates to the cen- 
tre). The star expands, moves back to the right in the HR 
diagram, and fills its Roche lobe after the C-burning con- 




time (yr) 

Figure 7. Case BB. Same as Fig. ^ for a 2.8 Mq helium star 
with a 1.4 Mq neutron star in a close orbit, during the case BB 
mass transfer. In the fifth panel, the upper and lower dotted lines 
show the CO and 0-Ne-Mg core masses (the latter is defined 
as the mass where the carbon abundance is below 0.1); and the 
shaded area mark the carbon-burning convective regions. 



vective core has disappeared (point b). During the second 
RLOF, Tcfi and L decrease while the period increases. After 
the first C-burning convective shell disappears, the mass- 
transfer rate reaches a maximum, Mmax,2 (point c), before 
decreasing again (see Fig. ^). 

In 3.1 - 5.8 Mq helium stars in intermediate orbits, 
RLOF takes place without detached stages. After Mmax.i, 
AI goes through a minimum and then increases again. Also 
due to the difference in thermal timescale between the core 
and the envelope, the occurence of the minimum AI (partic- 
ularly in this range of mass) varies. Like in the close-orbit 
systems, the period decreases throughout the whole mass- 
transfer phase. 



5.2.3 Wide-orbit systems 

Thermal-timescale RLOF is most rapid in wide-orbit sys- 
tems, and the maximum mass-loss rate, Mmax.i, is reached 
during the carbon core (or off-centre shell) burning. The 
mass-transfer rate then decreases and reaches a minimum 
value (in a 2.4 Mq helium star, the decrease in M leads to 
a detachment), before reaching another maximum. A/max, 2. 
The increases and decreases in AI can be seen in the HR 
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Figure 8. Case BB. Same as Fig. ^ for a 2.8 Mq helium star 
in an intermediate orbit. 

diagram as the up and down movements (see Fig. |^. Ex- 
cept in a 2.4 Mq helium star, RLOF in wide-orbit systems 
takes place without any detached stages. Again, only in 
Mhc < 2.9 Mq stars does the period reach a minimum 
value shortly after Afmax.i. 



pute the furthest (see the c notation in Table |^ in helium 
stars less massive than 4 Mq. Therefore, we expect that this 
is a typical phenomenon of the late stage of a case BB mass 
transfer from Mho < 4 Mq stars (in RLOF initiated during 
carbon shell burning (Habets 1986a), the envelope is thick 
so that the expansion of He-burning convective shell occurs 
deep inside the envelope). In Afnc > 4 Mq the start of a 
convective shell is only found in close orbits, but not in the 
wide orbits athough the calculations were done up to the 
same stage of evolution. However, we cannot exclude that 
such a convective shell does not eventually develop in these 
more massive stars as well. 

The outcome of mass transfer in the presence of a con- 
vective helium envelope is not clear. Mass transfer might be- 
come dynamically unstable which would lead to a CE and 
spiral-in phase. Although we do not find a convective enve- 
lope in helium stars more massive than 4 Mq in wider orbits, 
the large mass ratio [q ^ 2.6 at the onset of RLOF) and the 
high mass-loss rate (M > 10~* Mq yr"'^) suggest that a CE 
and spiral-in phase will very likely also take place. However, 
the high luminosity due to neutrino emission (L^ > lO'^ Lq) 
suggests that the core will collapse within a few decades af- 
ter the end of the calculation. The closeness to collapse is 
supported by the fact that Ne-burning has already started 
radiatively off-centre when the fourth C-burning convective 
shell appears. The explosion will probably be observed as a 
type Ic supernova. Nevertheless, should the two components 
undergo a spiral-in, a merger will not necessarily occur. A 
crude calculation of the gravitational binding energy of the 
convective envelope to the CO-core boundary of a 3.2 Mq 
helium star when the penetration takes place yields A = 
0.123. With J7CE = 1, Mho = 1.98 Mq^^ Mcorc = 1.62 Mq 
and Porb ~ 0'?22, we obtain from eq. (nol) a final separation 
of 0.0729 Rq (period of 1™89), which is still larger than the 
radius of the core (0.0275 Rq). Without helium left in the 
envelope, the explosion will be observed as a type Ic super- 
nova. The close post-SN separation subjects the two neutron 
stars to a strong loss of gravitational-wave radiation, caus- 
ing the two compact stars to merge on a timescale of 2 500 
years, which may produce a gamma-ray burst. 



5.2.^ The He-burning convective envelope and 
common-envelope evolution 

As can be seen in the right-hand panels of Fig. ^, a He- 
burning convective shell appears in the envelope before the 
appearance of the third C-burning convective shell. Such a 
shell also appears in previous helium star calculations (e.g. 
Habets 1986a, 1986b; Nomoto 1984, 1987; Pols 2002). Be- 
cause of the small envelope mass, this shell almost penetrates 
to the surface after the disappearance of the third C-burning 
convective shell. In response to mass loss from a convective 
envelope, the star expands, M increases enormously by 2 or- 
ders of magnitude, and the orbit widens. L and T^s decrease 
and the star moves to the lower right in the HR diagram (see 
Fig. ^). In a 3.8 Mq helium star, the He-burning convective 
shell appears and penetrates to the surface before the second 
C-burning convective shell appears. 

We are not always able to follow the evolution of the he- 
lium stars beyond the appearance of the second C-burning 
convective shell. However, we found that the behaviour oc- 
curs for those masses and orbits that we were able to com- 



5.2.5 Final masses and outcome 

The final masses and periods are shown in Table ^ Here 
final corresponds to the last stage of our calculations, which 
varies from the appearance of the first C-burning convective 
shell to the penetration of the He-burning convective shell to 
the surface. We found that after the carbon convective core 
has disappeared, only a small amount of mass is transferred 
to the neutron star. The period alters insignificantly. Hence, 
the values in Table |^ can be considered as the final masses 
and periods prior to SN ex plosio n, unless a spiral-in phase 
occurs as discussed in Sect. 



5.2.4 



With increasing initial period and decreasing initial 
mass, the amount of transferred matter decreases. In 2.4 
- 2.5 Mq helium stars, the mass of the remnant ranges 
from 1.3 - 1.6 Mq, but the mass of the CO core (1.21 - 
1.34 Mq) never exceeds the Chandrasekhar limiting mass, 
Mch- These stars will develop degenerate ONe cores. The 
future of these systems is not clear since they could not be 
computed through carbon burning. If the ONe core can grow 
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Figure 9. Case BB. Same as Fig. |^ for a 2.8 Mq helium star in a wide orbit. Instead of central abundances, the second panels show 
the stellar luminosity (dash-dotted-), and the contributions of helium burning (solid-), carbon burning (dotted-), and neutrino losses 
(dashed-line). The panels on the right-hand side show the evolution around the penetration of the He- burning convective shell to the 
surface. 



to 1.37 - 1.39 Mq, they will undergo a SN explosion trig- 
gered by electron capture on Ne and Mg (e.g. Nomoto 1984, 
1987), otherwise they will become heavy ONe white dwarfs. 

Stars with Mhc > 2.6 — 2.8 Mq which undergo case 
BB mass transfer develop CO cores more massive than Mch 
(where the actual lower mass limit probably depends on the 
orbital period). These stars will eventually complete neon, 
oxygen and silicon burning and form iron cores. The col- 
lapse of the Fe core by photodisintegration triggers a SN 
explosion, leaving a neutron-star remnant. 

The final amount of helium left in the envelope proba- 
bly determines whether the explosion will be a Type lb or 
a Type Ic supernova. The main observational criterion to 
distinguish between these types is the detection of helium 
in the spectrum. Helium is generally absent in SNe Ic, al- 
though a small amount may be present (Filippenko et al. 
1995). The maximum allowed amount of He is very uncer- 
tain, however, and may also depend on other factors such as 
mixing during the explosion (Woosley & Eastman 1997). 

For low- mass He stars and/or close orbits, we find there 
may be as little as 0.04 Mq of He left, and generally < 
0.10 Mq (see Table This is less than is achieved by 
stellar-wind mass loss in massive He-star models (Woosley, 



Langer & Weaver 1995; Wellstein & Langer 1999; Pols 2002) 
and also less than obtained from conservative case ABB or 
case BB mass transfer models by Wellstein & Langer (1999). 
We conclude that case BB mass transfer to a compact com- 
panion in a very-close binary is the most efficient way to pro- 
duce almost bare CO cores prior to explosion. These stars 
will almost certainly produce a Type Ic supernova. Such a 
model was first suggested by Nomoto et al. (1994) as the 
progenitor of the Type Ic SN 19941, and shown to match 
the observed lightcurve (Iwamoto et al. 1994). The amount 
of He left in the envelope increases with increasing initial 
mass and initial period. Probably, the more massive case 
BB remnants and those from the widest orbits undergo a 
Type lb explosion, as do helium stars that avoid case BB 
mass transfer (Pols 2002). 

5.3 The binary millisecond pulsars and double 
neutron-star pulsars 

The outcome of case BB mass transfer from 1.5 - 2.5 Mq 
helium stars are heavy white dwarfs with a neutron-star 
companion. Their positions in the period-mass plane (solid- 
and open-circles in Fig. cover the region where the exis- 
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Table 5. The orbital parameters of double neutron-star pulsars: 
total mass of the system in Mq, orbital period in days, and ec- 
centricity. 



PSR 


A/t 




e 


Reference 


B1913-I-16 
B1534-I-12 
J1518-I-49 


2.828 
2.679 
2.62 


0.323 
0.421 
8.634 


0.617 
0.274 
0.249 


Taylor & Weisberg 1989 

Wolszczan 1991 

Nice, Sayer & Taylor 1996 



tence of binary millisecond pulsars with heavy white-dwarf 
companions cannot be explained by extreme mass trans- 
fer (Tauris, van den Heuvel & Savonije 2000). Systems like 
J1756-5322 (Po,b = 0?453, Mwd = 0.683 Mq with the as- 
sumption that i — 60°) are not produced by our nor Tauris 
et al.'s calculations. Therefore, they must have come from 
stars with 5 < AfzAMS /Mq < 8. 

Case BB mass transfer from Mho > 2.6 — 2.8 Mq pro- 
duces neutron-star binaries. Attempts to find the progeni- 
tors of the three known galactic double neutron-star binaries 
(outside the globular cluster) have been carried out e.g. by 
Yamaoka, Shigeyama & Nomoto (1993), Fryer & Kalogera 
(1997), Bagot (1997), and recently FrancischeUi, Wijers & 
Brown (2001). Their calculations were approximated by the 
evolution of single helium stars, and it is often assumed that 
the helium star does not fill its Roche lobe (i.e. Rhc < Rl,Hc) 
in order to avoid an unstable phase of mass transfer. The 
orbital parameters of the observed binary pulsars are pre- 
sented in Table ^. Here we will try to see if our helium star 
calculations can explain their existence. 

If we simply assume that the radius of the pre-SN or- 
bit (the helium-star binary), flo, must be inside the radius 
turning point of the post-SN orbit, Op, (Flannery & van den 
Heuvel 1975) i.e.. 



(1 — e) ap < flo < (1 -f e) ( 



(11) 



where e is the eccentricity of the binary pulsar, then from 
Kepler's third law we have the following condition 



(1- 



- e)'M, 



P 

p 



< Mo Po < (1 + e) 



'Mp 



(12) 



where Mo, Mp are the total pre- and post-SN masses; Po, 
Pp are the pre- and post-SN periods. 

We plotted the allowed minimum and maximum pre- 
SN periods according to eq. ( |l2[ ) for each double neutron- 
star binary in Fig. llOl under the assumption that the pulsar 
companion is a 1.4 Mq neutron star. Fig. |l^ shows that 
systems like B1913-I-16 are easy to produce from systems 
with initial masses of 2.4 - 6.7 Mq in initial periods of 0'?1 - 
0^7, while systems like B1534-I-12 can be made from helium 
stars with the same range of mass in a narrower range of 
period (0'?3 - 0'?6). None of the case BB models we have 
computed can produce B1518-|-49-like systems. The remnant 
of the 2.4 Mq system (indicated by the arrow in Fig. px| ) is 
a ONe white dwarf /neutron star binary in a circular orbit. 
A case BC evolution might be able to explain the formation 
of this system, or a Mhg ^ 3 Mq star, which avoids RLOF 
altogether (the shaded area in Fig. [lo| ). 

The region to the left of the solid line in Fig. |l^ produces 
white dwarfs. The area below the solid curve will undergo 
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Figure 10. Range of mass and period of the immediate pre-SN 
progenitor of double neutron-star binaries; B1913-I-16 (dashed- 
), B1534-I-12 (dotted-), and J1518+49 (dash-dotted-line) . Upper 
and lower lines represent the maximum and minimum pre-SN pe- 
riods. The calculations done in this work are marked according 
to the type of remnant they produce: CO white dwarfs (solid- 
circle), ONe white dwarf (open-circle), or neutron stars (solid- 
stars). The orbital periods (Po) and masses Mhc f are those im- 
mediately prior to the SN or the formation of the white dwarf. 
The region below and to the left of the solid curve is the zone 
of avoidance for the formation of neutron stars. The shaded area 
marks the region where a double neutron star can be produced by 
avoiding RLOF, taken from Pols' (2002) single helium stars cal- 
culation after taking into account the effect of stellar wind mass 
loss. 




Figure 11. The final fate of the helium stars borm in systems 
consisting of helium star and neutron star: CO white dwarf (solid- 
circle), ONe white dwarf (open-circle), or neutron star (solid- 
star). Pi and Mho are the orbital period and helium star mass at 
the onset of the evolution of the binary. The lines separate the 
regions of case BA, BB, and BC mass transfers, taken from Pols' 
(2002) single helium stars calculations. The shaded area mark the 
region where mass transfer is dynamically unstable. 



case BA or dynamically unstable mass transfer (see Sect. y). 
Therefore, these regions can be considered as the zone of 
avoidance for the formation of double neutron stars. 
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6 CONCLUSIONS 

We have done calculations of helium stars with masses in 
the interval of 1.5 - 6.7 Mq with a 1.4 Mq neutron star 
companion. We allowed the helium star to fill its Roche lobe 
during helium core burning (case BA) and helium shell burn- 
ing (case BB mass transfer). The final fate of the helium star 
as a function of its initial mass and period is summarized in 
Fig. |ll|. 

Case BA mass transfer from helium stars with masses 
of 1.5 - 2.9 Mq occurs during helium core burning and stops 
when the central abundance of helium drops below 0.1 and 
the star contracts. In helium stars less massive than 2.4 Mq, 
RLOF is stable and takes place on the nuclear timescale. For 
2.4 < Mhc/Mq < 2.9, a rapid, thermal-timescale phase of 
RLOF is followed by a stable, slow phase of mass transfer 
on the nuclear timescale. Driven by radius expansion, a sta- 
ble phase of case BAB mass transfer occurs during helium 
shell burning. This phase of mass transfer stops when the 
degeneracy border is crossed and the star contracts. During 
RLOF, the entire envelope is removed from the helium star. 
The remnants of case BA mass transfer are heavy, degener- 
ate CO white dwarfs. This implies that a star with initial 
mass as large as 12 Mq still can become a white dwarf, as 
also found by e.g. Wellstein et al. (2001) in stars that go 
through conservative contact-free case A evolution. 

Stable case BB mass transfer from 1.5 - 2.1 Mq stars 
takes place on the nuclear timescale and stops when the 
degeneracy border is crossed and the star contracts. The 
entire envelope is removed during RLOF, and heavy, CO 
(from 1.5 - 1.8 Mq) or ONe (from 2.1 Mq helium stars) 
white dwarfs are the remnants in this range of mass. Case 
BB mass transfer from 2.4 - 2.5 Mq helium stars produces 
ONe white dwarfs and that from more massive ones pro- 
duces neutron stars. RLOF from Mne > 2.4 Mq takes place 
on the thermal timescale. In close-orbit systems, this is fol- 
lowed by a slow, stable phase of mass transfer. There is a 
tendency that the helium envelope becomes convective in 
helium stars less massive than 4 Mq. Mass transfer might 
become dynamically unstable in this case. However, the star 
probably collapses before a spiral-in can occur. Even if a CE 
and spiral-in phase does occur, the system is likely to sur- 
vive. In both cases, a double neutron-star binary is formed. 

From the amount of helium left in the envelope, we 
suggest that helium stars of high mass or in a wide orbit 
produce type lb SNe. Type Ic SNe come from helium stars 
of lower mass or in a close orbit. We are also able to provide a 
zone of avoidance for the formation of double neutron stars. 



ACKNOWLEDGEMENTS 

This work was sponsored by NWO Spinoza Grant 08-0 to 
E. P. J. van den Heuvel. It is a pleasure to thank Norbert 
Langer, Gijs Nelemans, Winardi Sutantyo, Thomas Tau- 
ris, and Lev Yungelson whose valuable comments and sug- 
gestions improve the early version of the manuscript; and 
Philipp Podsiadlowski for his constructive referee comments. 



REFERENCES 

Avila Reese V. A., 1993, Rev. Mex. Astron. Astrofis. 25, 79 



Bagot P., 1997, A&A 322, 533 

Bhattacharya D., van den Heuvel E. P. J., 1991, Phys. Rep. 203, 
1 

Brown G. E., 1995, ApJ 440, 270 

Chevalier R. A., 1993, ApJ 411, L33 

Chevalier R. A., 1996, ApJ 459, 322 

Delgado A. J., Thomas H.-C, 1981, A&A 96, 142 

Dewi J. D. M., Tauris T. M., 2000, A&A 360, 1043 

de Jager C, Nieuwenhuijzen H., van der Hucht K. A., 1988, A&AS 

72, 259 
de Kool M. 
Eggleton P. 
Eggleton P. 
Eggleton P. 
Eggelton P. 



351 
361 
279 



1990, ApJ 358, 189 
P., 1971, MNRAS 151, 
P., 1972, MNRAS 156, 
P., 1973, MNRAS 163, 
P., 1983, ApJ 268, 368 
Ergma E. V., Fedorova A. V., 1990, Ap&SS 163, 143 
Filippenko A. V. et al., 1995, ApJ 450, Lll 
Flannery B. P., van den Heuvel E. P. J., 1975, A&A 39, 61 
Francischelli G. J.. Wiiers R. A. M. J.. Brown G. E., 2001, 



sub- 



2001, A&A 376, 950 
1996, ApJ 466, L87 
A&A 231, 134 



mitted to ApJ, astro-ph/010321{ 
Fryer C, Kalogera V., 1997, ApJ 489, 244 
Habets G. M. H. J., 1986a, A&A 165, 95 
Habets G. M. H. J., 1986b, A&A 167, 61 

Hamann W.-R., Koesterke L., Wessolowski U., 1995, A&A 299, 
151 

Hamann W.-R., Schonberner D., Heber U., 1982, A&A 116, 273 
Iwamoto K., Nomoto K., Hofiich P., Yamaoka H., Kumagai S., 

Shigeyama T., 1994, ApJ 437, L115 
Landau L. D., Lifshltz E., 1958, The classical theory of fields, 

Pergamon Press, Oxford 
Nelemans G., van den Heuvel E. P. J. 
Nice D. J., Sayer R. W., Taylor J. H. 
Nieuwenhuijzen H., de Jager C, 1990, 
Nomoto K., 1984, ApJ 277, 791 
Nomoto K., 1987, ApJ 322, 206 

Nomoto K., Yamaoka H., Pols O. R., van den Heuvel E. P. J., 
Iwamoto K., Kumagai S., Shigeyama T., 1994, Nature 371, 
227 

Nugis T., Lamers H. J. G. L. M., 2000, A&A 360, 227 
Pols O. R. 1994, A&A 290, 119 
Pols O. R,. 2002, in preparation 

Pols O. R., Schroder K.-P., Hurley J. R., Tout C. A., Eggleton P. 

P. , 1998, MNRAS 298, 525 
Pols O. R., Tout C. A., Eggleton P. P., Han Z., 1995, MNRAS 

274, 964 

Savonije G. J., de Kool M., van den Heuvel E. P. J., 1986, A&A 
155, 51 

Soberman G. E., Phinney E. S., van den Heuvel E. P. J., 1997, 

A&A 327, 620 
Taam R. E., King A. R., Ritter H., 
Dewi J. D. M., 2001, 
Savonije G. J., 1999, 
van den Heuvel E. P. 



, 2000, ApJ 541, 
A&A 369, 170 
A&A 350, 928 
J., Savonije G. 



329 



J., 2000, ApJ 



Tauris T. M., 
Tauris T. M., 
Tauris T. M., 
530, L93 

Taylor J. H., Weisberg J. M., 1989, ApJ 345, 434 

Tutukov A. v., Fedorova A. V., 1990, Sov. Astron. 33, 606 

Tutukov A. v., Yungelson L., 1973, N. Info. 27, 70 

van den Heuvel E. P. J., 1994, in Shore S. N., Livio M., van den 

Heuvel E. P. J., eds.. Interacting Binaries, Springer, Berlin 
van Kerkwijk M. H., Charles P. A., Geballe T. R., et al., 1992, 

Nature 355, 703 

van Kerkwijk M. H., Geballe T. R., King D. L., van der Klis M., 

van Paradijs J., 1996, A&A 314, 521 
Webbink R. F., 1984, ApJ 277, 355 
Wellstein S., Langer N., 1999, A&A 350, 148 
Wellstein S., Langer N., Braun H., 2001, A&A 369, 939 
Wolszczan A., 1991, Nature 350, 688 

Woosley S. E., Eastman R. G., 1997, in Ruiz-Lapuente P., Canal 



© RAS, MNRAS 000, 0-|l| 



16 Dewi et al. 

R., Isern J., eds., Thermonuclear Supernovae, NATO ASIC 

Proc. 486, p. 821 
Woosley S. E., Langer N., Weaver T. A., 1995, ApJ 448, 315 
Yamaoka H., Shigeyama T., Nomoto K., 1993, A&A 267, 433 



